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^ ' Abstract 

^ ; 

^ We consider variations of coupling strengths and mass ratios in and beyond 
the Standard Model, in the light of various mechanisms of mass genera- 

O^I I tion. In four-dimensional unified models, heavy quark and superparticle 

■ thresholds and the electron mass can completely alter the (testable) rela- 
I tion between A In a and Aln/i, where n = mp/rrie- In extra-dimensional 

■ models where a compactification scale below the fundamental scale is vary- 
l/^ \ ing, definite predictions may result even without unification; we examine 

■ some models with Scherk-Schwarz supersymmetry-breaking. 

m ■ 
O ; 

^ ■ 1 Introduction 

Q^' If the recently measured cosmological variation in a T is to be explained within 

1^ ■ a unified model, i. e. one in which all couplings are constrained by a relation 

. . \ at a certain energy scale, definite predictions should ensue for the variations 

' of other quantities accessible to astrophysical measurement, in particular /x = 

^ ■ mp/rrif. 3 and gp (the gyromagnetic ratio of the proton) [^[^J. The "theory 

of everything" which would be required to make such predictions does not at 
present exist, although there are candidates for large parts of it. So far, results 
have been somewhat discouraging for unification, even considering the amount 
of theoretical uncertainty attached ^21 El |H1 El IH^ ■ 

To summarize the situation very roughly, the mechanism(s) of mass genera- 
tion in d = 4 field theory, including dimensional transmutation in QCD as well 
as various ways of generating a Higgs v. e. v. and fermion mass hierarchies, may 
be very sensitive to changing couplings at high scale, since they can involve 
slow logarithmic running into the strong coupling regime. But the running 
of the U(l) and SU(2) gauge couplings is weak and, crucially, a does not de- 
pend exponentially on a gauge coupling at high scale, but only linearly (to first 
approximation) . 

Thus the fractional variations in the mass ratios which determine // (and may 
affect gp) obtain contributions many times the size of the fractional variation 
of a. However experimental bounds on variations in /i and gp are of a similar 
order to the signal in a, namely a few times 10~^. Hence the number of theories 
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that could possibly be consistent with experiment appears drastically reduced 
and may indicate a fine-tuning problem. Such a conclusion turns out to hold 
also if the parameters of the Higgs and flavour sectors are held fixed (relative 
to the fundamental scale), since the variation of the SU(3) coupling at high 
scale is almost always at least as large as that of the U(l) and SU(2) couplings 
(whether the variation arises from a varying unified coupling constant, or a 
varying GUT scale relative to the fundamental scale, or some combination) 
and the exponential amplification due to strong running of QCD is always 
present [El [El CEl EH • 

Our aim is not to cure this apparent problem, since it is obviously desirable 
that theories be ruled out by experiment. Instead we will take more fully into 
account effects associated with thresholds in RG running and mass generation, 
and show that they may significantly alter theoretical expectations, and also 
point out that predictions can be made in the case of a varying compactification 
radius even in models without unification. 

Important thresholds arise already in the SM, for example from the 
and heavy quarks; theories capable of addressing the hierarchy problem are 
expected to possess additional thresholds, and imply particular relations for 
variations in the mass spectrum (which, of course, directly influences //). To 
make predictions, one must identify a single varying degree of freedom in the 
model as the source of the observed effect: given the complexity of realistic 
unified theories there may be more than one candidate for the source of the 
variation. Nevertheless certain possibilities are obviously simpler and better- 
motivated than others. 

For a varying GUT gauge coupling, given the SUSY-GUT value Mq — 
2 X 10^6 GeV and without considering variation of mass ratios in the electroweak 
and flavour sectors ^lEllBj, earlier calculations gave R = Aln/i/Alna = 
36 with an uncertainty of a few due to higher-order effects. Our estimates, 
including varying mass ratios ^ , and with an improved treatment of the proton 
mass ^ give ^ = —13 it 7 in a theory with softly-broken supersymmetry at the 
weak scale and = 4 it 5 without. 

1.1 New variations and twists from extra dimensions 

Various types of model with extra dimensions |321 1331 14U1 1411 1431 1441 145j may 

have interestingly different predictions since first, gauge and matter fields may 
propagate in different submanifolds in extra dimensions; second, the RG scale 
at which the variation of couplings originates is not tied to any unification scale; 
and lastly, the mechanisms of breaking SU(2)xU(l), supersymmetry or GUT 
groups in some theories, have a distinctive dependence on the radii of compact 
extra dimensions. 

On the first point, traditional Kaluza-Klein-like theories have a universal 
dependence on the size of the extra dimension, since all fields live on the same 
manifold; however, "brane-world" theories restrict some fields to a submanifold. 

'^In the case when electroweak symmetry-breaking occurs in a hidden sector with a gauge 
couphng varying with the unified coupling 
^See also 
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The submanifold may indeed only be x S^, in which case couphngs involving 
localized fields are not affected when the size of transverse dimensions varies. Of 
course, if all fields live on the same brane, there is no new behaviour of varying 
couplings (except for gravity) resulting from the presence of extra dimensions. 

On the second point, the couplings need no longer be determined (by, say, 
a dilaton v. e. v.) at GUT or string energy scales and run through 14 or more 
orders of magnitude to reach experimentally-accessible energy scales. Both the 
fundamental cutoff scale, and compactification scales, can be much lower, thus 
RG running may be quantitatively different: the scale at which the variation 
in couplings originates is different. This possibility includes the case of extra- 
dimensional GUTs with Kaluza-Klein thresholds |33| I34j . In fact, even non- 
GUT theories may have well-defined predictions, in the case where the variation 
is due to a cosmologically varying radius. We only need to know the physics at 
and below the energy scale where the variation is transmitted to the SM fields: 
in general couplings at this scale will not be unified. (Compare |21j . where the 
unification relations are irrelevant for the effect discussed). 

On the third point, Scherk-Schwarz (S-S) symmetry breaking [3511321 HOI I41j 
generates a mass scale inversely proportional to the radius of compactification, 
while perturbative couplings in the effective D = 4 theory vary as a power of 
the radius (depending on which fields propagate in extra dimensions ^). Thus 
for a varying radius, mass ratios and coupling strengths (defined at the scale 
of the inverse radius) will both vary as power laws with small exponents. This 
contrasts with the exponential dependence of mass scales on the high-energy 
gauge coupling in the case of mass generation by dimensional transmutation in 

1.2 Other issues 

There are many aspects of varying couplings that we will not discuss in detail. 
One outstanding question is what form the time- (or space-) dependence takes 
and what dynamics generate it. Following on from this question, if the variation 
originates from a light scalar, long-range forces and apparent violations of the 
weak equivalence principle will likely ensue HSl EUl UHl '■ see also Section 13.21 

If the spacetime dependence could be predicted, then bounds from the Oklo 
phenomenon [Bi and from the early Universe (in particular nucleosynthesis 
and CMB) [Sl[3|22] might well be the most stringent.^ We will take the opposite 
approach and concentrate on the range of redshifts at which a nonzero variation 
in alpha is claimed: direct constraints arise from astrophysical observations at 
the same epoch. See ^1 El El jSl] for similar approaches. 

^Classically, the 4d coupling constant of a gauge field q-*' propagating in extra dimensions 
varies simply as , where Ri is the radius associated with a one-cycle of a compactification 
manifold, which in the simplest case has the metric of a torus. 

*One cannot alter the exponential dependence of Ac, the confinement scale of strong inter- 
actions, on the asymptotic SU(3) coupling constant without radically altering phenomenology. 

^These bounds, which depend on nuclear or gravitational physics, do not directly limit 
a, rather they probe some combination of a and other dimensionless couplings: thus one 
can imagine degeneracies, enhancements or even cancellations depending on how other 

couplings may vary. 
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The case of spatial variation [HJ is less straightforward. Binning points by 
redshift is then not useful, and the dependence of alpha on the source (actually, 
on the absorbing medium) could be highly model-dependent. It could be that 
couplings measured on Earth are different from those in the gas clouds which are 
under observation simply by virtue of the different astrophysical environment 
- possibly due to the different gravitational potentials OUj . 

Whatever the dynamics producing the variation, one can investigate corre- 
lations between (variations in) the quantities that determine the spectral lines, 
a, Qp and /i. That these quantities are not measured simultaneously for the 
same sources or absorbers may be of concern. In models with variation over 
cosmological time, we have to assume that the variation is smooth and rela- 
tively slow. In models with spatial variation, we might expect a correlation 
between the variation of a in one source or absorber, relative to on Earth, and 
variations of gp and /x in another; however in the case of strongly environment- 
dependent variation we would need much more data on the latter quantities, 
which are currently bounded by a small handful of systems. If the spatial 
variation is long-range, rather than environment-dependent, one can say little 
without more detailed work on the spatial distribution of the sources. Spatial 
variation might be diagnosed by measurements at low redshift, or by analysis 
of spatial correlations (requiring a much larger data set than at present). 

The problem of fine-tuning the mass of any scalar field responsible for the 
variation against radiative corrections is an acute one, as noted in TT (in 
the context of quintessence) and ^S|- we have no further insights into this. 
The possibility noted in |2J of a rather heavier scalar m ~ 10"'^ eV in a new 
weakly-coupled sector alleviates such a problem, but the gravitational energy 
of the scalar must be dynamically cancelled if the model is to be cosmologically 
viable. 

Finally, the possibility of degeneracy in the interpretation of measurements 
in terms of atomic physics has been raised Jl]: the changes in spectral lines 
interpreted as varying alpha could in principle be due to varying /i or gp instead. 
A full evaluation for the many-multiplet technique could be complicated, as one 
needs to know the dependence of many spectral lines on all three parameters; 
recall that numerical calculations were required to find the expected variations 
in frequency due to alpha alone. We will, for the moment, take the data on a, 
fj, and gp at face value. 

1.3 Note added 

After the first version of this paper appeared. Murphy et al. released a preprint 
^ with a revised value of Aa/a, with the change being due to the identification 
of an additional source of random scatter in some of the systems studied, and a 
corresponding increase in the statistical errors assigned. This new value is used 
in the discussion in Section^ 
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2 Mass generation and thresholds 



The importance of the Higgs v. e. v. (or the ratio vn/Mp) was realised in [5], 
which found that its dependence on the SUSY-breaking masses and the top 
Yukawa was dominant in determining the outcome of nucleosynthesis. We will 
find that the variation of vh (and superpartner masses) can also be decisive 
for the spectral data, since they depend crucially on mass ratios, particularly 
the ratios of quark and lepton masses to the scale of QCD. After a general 
discussion of how thresholds enter the low-energy observables, we give examples 
of how the inclusion of thresholds, formally a higher- loop effect, can completely 
alter predictions from unified theories. 

We fix notation by quoting the (solution of the) one-loop RG equation for 
gauge couplings: 

«rH^-) = «rH/^+)-l^infeV (1) 



2tt 

thus Pi is negative for asymptotically free groups. Now given a charged field 
whose decoupling mass is m*'^(m*'^) (in a convenient mass- independent renor- 
malisation scheme), we have 

„-.<,-,.„,H.^,-|.(^)-|,„(^) 

where bf^ = bf — b~ , the beta- function coefficient being bf above the thresh- 
old and b~ below, with tree-level matching at m^^. In the case of multiple 
thresholds one sums the corrections to ■ For the QCD invariant scale 
Ac = Me~^'^/^"^^*^^ where < M < rric we find for non-super symmetric 
theories 



- (3) 



where /x"*" > mt, and for superpartners (squarks of geometric average mass niq 
and gluinos of mass rrig) 



2/27 / \ 2/9 



(4) 



where ix^ > mgrrig and the powers 2/27 and 2/9 come from ratios of beta- 
function coefficients, specifically —b^^/b^^"^"^ = ^3^/9. Thus the variation in Ac 
is 

/^"^^Ac ^ 2tt Aa3(^+) ^ 2 + ^ /^/i^A^ + — 

Ac /x+ 9a3(/j,+) a-i{ii+) 27 ^^^rriq fi+ 9 \mg /i+ nig fi+ 

(5) 

where terms in brackets are to be ignored for the non-supersymmetric case. 
The heavy quark masses to be used are strictly the decoupling masses mq{mq), 
but we showed in that the difference between using this definition and the 
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"invariant quark mass" rhq defined such that it has no scale dependence, is 
neghgible. 

We see immediately that the threshold terms are of higher order in 03 , thus 
formally they ought to be grouped with the power-law correction to Ac from 
two-loop running. However, when the Higgs v. e. v. is varying rapidly compared 
to a, as is generic in theories unified at high scale jl5l 15]. the threshold terms 
can become dominant, in contrast to the two-loop term, which is a model- 
independent infrared effect. 



2.1 Electroweak symmetry-breaking 

One case where thresholds are important is when electroweak symmetry-breaking 
(EWSB) is triggered by nonperturbative gauge theory effects and the gauge 
group giving rise to such effects is unified with the SM gauge couplings at high 
scale. This case encompasses gravity- and gauge-mediated SUSY-breaking and 
technicolor or composite theories, under the assumptions that all gauge groups 
be unified, and all mass scales lower than the fundamental scale Mx (which may 
be the GUT scale, the string scale related to a' or the Planck scale of whatever 
underlying theory we are considering) are generated dynamically. The generic 
dependence of vh is 

^ = A:a^e-2™/'"^"^ (6) 
Mx 

where fc is a numerical constant and n parameterises a possible power-law de- 
pendence, for example through the Higgs quartic coupling which may (as in the 
MSSM) be related to a gauge coupling varying with ax, or through the gravi- 
tational constant in the case of gravity-mediated SUSY-breaking. The value of 
m is model-dependent, however the product m/bhOix can be estimated more re- 
liably. In the case of theories without fundamental Higgs, this equation simply 
parameterises the v. e. v. of whatever condensate breaks SU(2)xU(l). Thus, 

Mx A Vh ( , 27rm ^ Aax 

ax 



+ -. (7) 



Now under the assumption that k and n are order 1, we neglect (the logarithm 
of) their contribution to Eq. © and find 

~ \n[Mx/vH ~ (2 X 10^'^)/(2 x 10^)] ~ 32 

hhax 

where we identified Mx with the SUSY-GUT value. ^ Thus in Eq. ® the RHS 
becomes approximately {n + 32)Aax/ax; which is consistent with neglecting 
n of order 1. If instead we take the heterotic string scale Mx ~ 4 x 10^^ GeV 
we evaluate the RHS as (n -|- 35), thus the uncertainty introduced by a small 
integer value of n can be neglected. For illustrative purposes we will take 

®In the non-supersymmetric case, single-step unification is of course inconsistent witli data. 
Witli tfie nominal values Mq — 10^^ GeV and ~ 42 we obtain 29 rather than 32 for 
this equation. Non-SUSY unification is of course possible with intermediate breaking scales, 
however such theories lose predictivity. In what follows the "non-supersymmetric" case will 
just mean the result when superpartner thresholds are ignored. 
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lS.hi{vH /Mx) = (34 lb 2)Alnax- From standard GUT relations, neglecting 
the effect of thresholds on the running of SU(2) and U(l)y one finds ^1 E] 
Alna = Alnaem('^e) — O.SAlnax, thus we may also write /\hi{vH /Mx) = 
(68 ± 4) A In a J Similarly, we may parameterise SUSY-breaking masses as 



fc'a^e-^^"'/*''"^ (8) 
Mx 

where rn-g, nig etc. vary as m, thus their variation can also be estimated as 
Aln(m/Mx) = (68 it 4)Alna given that superpartners are around the elec- 
troweak scale. 

It has been pointed out ^1 IH] that radiative EWSB may depend sensitively 
on a combination of the top Yukawa coupling and ratios of SUSY-breaking 
soft mass terms. Hence, any relation between vh/Mx and fh/Mx may be 
significantly more complicated. Variations in yt and soft term ratios are ex- 
ceedingly model-dependent, to the extent that it is futile to look at detailed 
models. The best we can do is to parameterise our ignorance (somewhat 
as in ,15^) by allowing the Higgs v. e. v. to have a different fractional varia- 
tion Aln(fj|f/Mx)/Alnax compared to that of the average superpartner mass 
Eq. We put 

_^\n{vH/Mx) . _Aln(m/Mx) 
Pv = T-; ' PS = 



A In ax A In ax 

where (3^ is defined for all models whether SUSY or non-SUSY and (3$ only 
in the supersymmetric case. Thus, from the above discussion, if there is high- 
scale unification and EWSB is generated by dimensional transmutation, we 
expect j3y ~ \i\{Mx /vh)-, which also holds for the SUSY case in the absence of 
strong dependence on yt and soft mass ratios; if SUSY-breaking is triggered by 
nonperturbative effects in a hidden sector we have f3s — ln(Mx/?Ti). 

If the hidden sector group is not unified with the GUT group, in the sense of 
having the same coupling constant at Mx, these expectations may be altered. 
In some heterotic string models one obtains gauge kinetic functions S + eT 
for the visible sector and S — eT for the hidden sector, where S is the (four- 
dimensional) dilaton and T a volume modulus, and 1 > e > 0. If we identify 
S as the varying d. o. f. the model behaves similarly to the case of a "unified 
hidden sector" ; if we identify T then the values of j3s and (3y will be large and 
negative! Hence, we keep (3$ and (3^ as free parameters until the final stage of 
calculation. 

2.2 Initial estimate of /i 

In the first approximation we neglect variations in Yukawa couplings at high 
scale, and obtain 

Aln^ - /^AAlna3(M.)=(^^ + ^/?. + l/3.)Alna3(M.) 

= (24 ± few) [non-SUSY], (39 ± few) [SUSY] (9) 

Ox Oix 



^See section 031 for a treatment including thresholds. 
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using 1/ax — 24, where we estimate "few" as about 3 and set /J^ = for a 
theory without superpartners. If the proton mass is weh-approximated by a 
constant times A^, we can estimate 



^ = (-10 ± few)^^ [non-SUSY], (5 ± few)^^ [SUSY]. (10) 

H ax ctx 

The dependence of the electron Yukawa couphng may also come into the esti- 
mate, possibly at the level of changing this value to, say, —13 [+2] if we imag- 
ine Uf, oc a\ due to some dynamics of flavour structure. This expectation 
contrasts with the cases where the Higgs v. e. v. and superpartner masses are 
fixed relative to Mx, giving A/i//i ~ 17 Aax / oix, or where vh varies as in 
(O but the resulting effects on and of QCD thresholds are neglected, giving 
A/x//x ~ —n/S.ax/oix- Thus the mechanisms of mass generation and the pres- 
ence of thresholds, especially superpartner thresholds, can completely alter ex- 
pectations for the relations between varying // and ax- In the next section we 
look at the observables a and /x in more detail and obtain preciser estimates 
taking mass generation and threshold effects fully into account. 

2.3 More detailed treatment of a and mp 
2.3.1 a and thresholds 

At first sight, "varying the unification scale" and "varying the electron mass" 
while keeping other parameters "constant" should produce equivalent effects 
on a: only the ratio of two masses is physical. However, the behaviour of 
intermediate thresholds is most important to determine the result. Hence, we 
include all known charged thresholds to determine what effect they may have 
on naive GUT predictions. 

"Varying Mq^ is interpreted as all low-energy threshold masses (that orig- 
inate from EWSB) varying by the same fraction with respect to Mq'- thus 
m^ji/MG, nir/MG, my/ /Mq, etc. vary by the same fraction as m^jMG. ^ "Vary- 
ing rrie is interpreted as holding all threshold masses fixed in units of Mq except 
for rUe/MQ itself. The beta-function coefficient per Dirac fermion of charge Q 
IS 4QV3, thus "varying me" leads to the relation 

a-i = aT\^+) - In , _ = A In ^ « 7.7 x lO'^ A In 

where mg < /i^ and we take a^^{fj,~^) to be fixed. 

Since we have a~^ = ai^+a2^ at the weak scale and TrQ^ = Tr Tg +y^Tr 1 
over an SU (2) irrep, we easily see that the effect on a of a certain fractional 
variation in the masses of a SU (2) multiplet with mass below Mw is no different 
from the same fractional variation in a multiplet lying above Myy. We have in 
general 

Aa _ ^ a Aai(/i+) ^ Q^^^ m^^ 



« it!^2"*(^^) ai(^+) ^ 27r /i^ 



+ „V^^-^Aln^ (11) 



Since the GUT scale is defined through the masses of superheavy fields that trigger the 
breaking to Gsm, their threshold masses are constants in GUT units. However, see |16| and 
the discussion preceding Section 3. 
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where ai = g''^ /Att, the second sum is over all charged fields, f^^ is 2/3 per 
chiral (or Majorana) fermion, 1/3 per complex scalar and 11/3 per vector bo- 
son. All logarithms ln(m*'^//U^) are numerically comparable if fx^ is very large, 
however in the case of low-scale models the difference between ln(me/Mx) and 
\x].{myi/ /Mx) may be significant. For the light quarks, the confinement scale 
provides a dynamical cutoff. We have, setting aside the first term in Eq. (|11() . 

^Ot 1 f ^ o \ ^ ^ f ^ \ ^ l^cfnt 1 ^ 1 TTlb ^ ^ , "^1 

|th 



itn - , -^AAlnax + - -Aln „ -F-Aln — ^-h3Aln 

a' 137-27rV3 3 Vs Mj. 2, Mx Mx 

—A In — ^ -F8Aln h-Aln— ^ (12 

1 /8 2^ \ 

r/?A + /3. + — /35 Alnox (13) 



137 - 27r V3 3 

where mi is an averaged lepton mass, we include the SM Higgs contribution 
as the longitudinal modes of the VF's and count the remaining modes of the 
second Higgs doublet in SUSY models as if superpartners. m and rriH stand for 
geometric averages over superpartners and heavy Higgses respectively, where 
such fields exist. 

Now in the case where ax varies at a high scale, taking /3„ ~ = 34 it 2 
and using the dependence of Ac from Eq. © both with and without SUSY 
thresholds we find 

— |th=i (0.11 ±0.01)^^ [non-SUSY], (0.49 ± 0.03)^^ [SUSY] 

a ax otx 

which is a non-negligible correction to the direct contribution from varying 
ai,2(Mx), (Aa/a)direct = (8a/3ax)Aax/ax — 0.47Aax/ax- (With the the 
non-SUSY SU(5) estimate ax ^ 1/42, this ratio is 0.82.) 

2.3.2 nip contributions from light quarks and electromagnetism 

The proton mass receives several contributions: studies in chiral perturbation 
theory |3H| indicate a dominant contribution from nonperturbative effects in 
the chiral limit, which can only be proportional to Ac, and subleading terms 
depending on n, d, and s quark masses. Heavy quarks are assumed to be well- 
described by their threshold effects on Ac. There is also an electromagnetic 
contribution 132] • 

The effect of varying light quark masses can be found by applying the 
Feynman-Hellmann theorem: 

ruq nip mq dnip / , - , \ 

— A— ^ ~ = niq{p\qq\p) 

vrip niq nip oniq 

where the matrix element is essentially the so-called cj-term. Quark masses 
appear in this equation such that any multiplicative change of normalization of 
the q operators cancels out: physical quantities do not depend on a particular 
definition of quark masses. Hence, as for the heavy quarks, we use asymptotic 
or "invariant" quark masses directly proportional to the Higgs v. e. v. times a 
Yukawa coupling. 
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In the limit of isospin symmetry we have 

A In ^ = m^V^7v(0) = 0.048 ± 0.01 

m ^ 

where m = {mu+md)/2, using CTtt at (0) = 45±10MeV j^EH]- Isospin-violating 
effects can be shown to be small by considering pion-nucleon scattering |^ or 
nip — rUn, hence we use the coefficient 0.048 for both u and d quarks. 

The strange contribution is related to the strangeness content of the nucleon 
y as 

rris drup _ msy anNjO) ^ 2(p|ss|p) 
mpdnis m 2 nip ' {p\uu + dd\p) 

An estimate dlnnip/dlnnis ~ 0.2 was made in |22| based on the lattice result 
{p\ss\p) ~ 1.5 with a sea-quark mass of 0.154MeV [SJ, corresponding to y = 
0.36 lb 0.03. However this large value was contested in |52j where negative y 
consistent with zero was obtained (negative values are unphysical). A safe but 
imprecise estimate is thus to take d In nip /d In nis = 0.12 it 0.12, corresponding 
toy = 0.2 ±0.2 [IHl. 

The electromagnetic contribution, composed of the direct electromagnetic 
self-energy and a part originating from nondegeneracy of pion masses |49| , gives 
nip^dnip/da ~ 1.8 x 10~^. Except in the case of nucleosynthesis, isospin- 
violating and electromagnetic effects can be neglected.^ For consistency we 
must put 

A dnip_^ y ^ 0.78 i 1 



m„ 9A dlnniq 

q=u,d,s ^ 



Putting it all together, we have 

A , I 2-K Aa3(M) 2 \ - . , TTig 4 . , m 
Aln^ ~ 0.78 — +-T^H > Aln— ^ +-Aln — 

M I 9a3(M) asiM) 27^^ M 9 M 

+ (0.12 ±0.1) A In ^+ 0.048 V A ln^ + --- 

u,d 

^•^^ -A In a3(M) ± 0.058 V A In ^ ± 0.35A In — 



a3(M) ^ ^ 

± (0.12 ±0.1) A In ^± 0.048 V A ln^±--- (14) 
M ^ M 

u,d 

where apart from the strange term the individual coefficients have errors of 
order 10%. Now as before we can take the simplest unification scenario in 
which all superpartner and quark masses vary approximately with vh/Mx to 
obtain 

777 /O 54 \ 

Aln^ ~ ^ + (0.39±0.12)/3„ ±0.35/3s Alnax (15) 
M \ax J 

~ (13± (13±4) ± (12± 1)) Alnax 



®For a consistent unit-free treatment of nin — mp in nucleosynthesis allowing both a and 
QCD parameters to vary, see |15|: see also |30| on implications for composition-dependent 
forces. 
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where the last term on the R.H.S. is to be discarded in the absence of super- 
partners. Thus the inclusion of light quarks and SUSY thresholds and con- 
sequent reduction in the contribution has a non- negligible effect. Using 
Aln/i/Alnax — (Aln(mp/Mx) — /3^AlnQ;x)/Alnax we come to the result 



Without variation in the weak and superpartner scales (setting (3v = (3s = 0) 
we arrive at the rather model-independent expectation ^ = 25 it (few). This 
differs somewhat from the expectation ^ ~ 36 based on rup = const, x Ac and 
ignoring the effect on a of varying Ac through light quark thresholds. (However 
we will see that this scenario is still ruled out.) 

If we allow the weak scale to vary but ignore SUSY thresholds, setting 
I3y ~ 34, ax = 1/24 and Ps = 0, we obtain ^ = — 13 it 7. The change here 
is mainly due to the variation in rrie/Mx- If we now change (5$ to 34 it 2 we 
arrive at ^ = 4 it 5. Here, Aln/i turns positive and A In a becomes somewhat 
larger through the variation of SUSY thresholds. 

The "varying Mc" scenario of for which ac is constant can be recovered 
by taking the limit /S^, — > oo. In the presence of one should send (3$ to infinity 
and also enforce (3.^ = f3s- The expectation in this case is for R = —330 ± 65 
without superpartner contributions, or R = —20 it 9 with. Again, these differ 
from previous estimates (13) due to our more detailed treatment of mp and of 
light quark thresholds. The enormous value of R in the "non-SUSY" case is 
due to a cancellation in Aa in the absence of the "direct" contribution. 

The string-inspired scenario of is more complex since it posits a string 
scale Mst at which the gauge coupling is invariant, as well as an effective GUT 
scale Mg varying relative to Mst- SUSY-breaking and EWSB are taken to be 
independent of Mc/Mst- There are two ways to implement the scenario: firstly, 
one may identify Mx with M^t and set (3^ = [3s = ^ = ctx] however one has 
to rederive Eqs. (jHHSI)! (fT2HT3|l . etc., explicitly including heavy GUT thresholds 
whose contributions oc ft-'^A ln(MG'/Mst) can be found given the beta-function 
of the unified group. 

The second way is to identify Mx = Mq and then impose f3v = Ps = 
— Aln(MG/Mst)/AlnaG' (since the weak and supersymmetry scales are invari- 
ant relative to Mgt). The variation in ac is found by RG running down from 
Mst to Mq. The results would differ slightly from those of again due to 
different treatment of rUp and light quark effects, but the overall conclusion 
that (even in the absence of rapidly-varying weak and SUSY scales) the great 
majority of possible GUTs result in values |^| > 10, should remain. 

3 Varying radii of extra dimensions 
3.1 Rationale of a varying radius 

In a model which is described by (4 + (5)-dimensional field theory over some 
range of energies, there is a hierarchy between the (inverse) compactification 



A In a 



Aln^u 



= R 



OMa^ + (-0.61 ± 0.12)/?^ + 0.35^5 
0m2aZ^ + 0.0018/3^ + 0.011/?s 



(16) 
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radius or radii R~ and the ultraviolet cutoff of the higher-dimensional theory 
A/), where D = 4: + 6. This makes it a well-defined problem to compute the 
effect of varying radius on the 4D low-energy theory, since one can take physics 
above the scale R^^ to be unchanged by such variation (apart from the change 
in masses of Kaluza-Klein modes) , consistent with decoupling. It is conceivable 
that physics above this energy scale is also varying, but such variation would 
likely not be calculable since it might involve physics at energies equal to or 
above the cutoff. 

The superstring dilaton and moduli are a priori candidates for cosmologi- 
cally varying scalar fields and may be light enough to survive at low energies, 
yet they parameterize inverse radii which may be in the range 10^^~^^ GeV. In 
principle they might come into our analysis, except that the theory above the 
inverse radius scale may be "stringy" to a significant extent: the cutoff is at 
the string scale (or near the 11-dimensional Planck scale for M-theory), hence 
if the inverse radius also approaches this scale there is no range of energy over 
which extra-dimensional field theory makes sense. In perturbative string theory 
the dilaton- and moduli-dependence of low-energy quantities are extractable, 
but we face the usual problem that string models which are realistic enough to 
contain a reasonable approximation to the SM spectrum, are too complicated 
for detailed calculations to be practicable. One can however make simplifying 
assumptions based on the behaviour of less realistic models, leading to results 
such as those discussed in [HI IH1 • 

3.2 Are long-range forces inevitable? 

Naively, one might expect that a varying radius which was the source of varying 
alpha would mediate observable long-range forces. The reasoning is that, given 
that the scalar field appearing in the four-dimensional effective action ("radial 
modulus" ) is evolving in a potential on timescales of order the inverse Hubble 
rate, its mass must be comparable to the Hubble rate, thus it mediates a force 
over very large distances. The strength of the force (i. e. its coupling to matter) 
in Kaluza-Klein type theories is found at tree level to be comparable to gravity 
(see e.g. 24 ) thus it is likely ruled out by Solar System tests of deviations 
from general relativity and tests of the Weak Equivalence Principle ^2F. (Of 
course one would not observe any deviation from the inverse square law in the 
laboratory.) 

However, this reasoning is based on a number of assumptions which are not 
universally valid: therefore although it appears difficult for a varying radius to 
evade these bounds, it cannot be dismissed immediately. 

The first assumption is that the apparent measured variation is due to a 
smooth evolution over billions of years. However, as we mentioned in Section 
11.21 the variation may be spatial and long-range or short-range; the variation 
may even operate on distance scales as small as a galactic halo. Hence the 
mass of a varying scalar may be some orders of magnitude above the inverse 

^'^Even if one is allowed to adjust by hand the scalar kinetic term and couplings to matter, 
if one specifies that it should be the origin of the alpha variation, then apparent violations of 
the WEP should be not far below current bounds |28ii30( . 
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Hubble radius. This does not much help to evade observational bounds, but it 
is worthwhile to refute an unwarranted assumption. 

The second is that the radial modulus is identified with the light scalar 
that is slowly rolling in a non-varying potential. However, compactified radii 
may vary without the modulus being light. There is an analogous argument 
for two scalar degrees of freedom in the SM, the sigma meson corresponding 
to excitations of the chiral condensate and the Higgs. In the previous parts 
of this paper we saw that these scalar v. e. v.'s are likely varying (in Planck 
units) in a way correlated with alpha. By the same argument one would also 
conclude that the sigma and the Higgs should have masses of order the inverse 
Hubble length. In fact, the scalars are massive and always at the minima of 
their potentials: their variation is induced by a varying potential. There is 
no reason why a variation of the radius should not similarly be induced by 
a varying stabilization potential, with the mass of the modulus being 10-3 eV 
(the current limit of sensitivity of fifth-force experiments) or greater. We do not 
specify the source of such a varying potential, except that it should not directly 
affect the coupling constants of SM fields. Thus the radius may mediate the 
variation to the observable sector. 

The third is that the scalar couplings to matter are constant. However, even 
in the simplest case of unification, this is not the case, since any variation in the 
SU(3) coupling at high scale results in nonperturbative dynamics at low energy 
which determines the masses of nucleons as an exponential function. When one 
combines this fact with the presence of Scherk-Schwarz SUSY-breaking, elec- 
troweak symmetry-breaking, varying Yukawa couplings and other low-energy 
phenomena it becomes clear that the effective couplings to matter will be field- 
dependent, possibly in a complicated way. It was shown in j^H] that, in the case 
of field-dependent couplings to matter, if there is a value of the scalar where 
the coupling is small, this value is an attractor for the evolution of the scalar. 
Hence the deviations from GR may be suppressed. This mechanism was dis- 
cussed in a model with varying radius in |27] , although the details of this model 
are different from those I discuss below. 

The fourth is that the tree-level action for the radial modulus is a good 
approximation. However, it is possible that there are large corrections to the 
kinetic term and couplings. In |27j radiative corrections to the action were 
mentioned as a possible way to evade bounds on a light scalar. 

Violations of the WEP and other deviations from GR are very strong bounds 
on all models which seek to explain the dynamics of varying alpha, not only 
models with varying radii. We do not have any definite solution to the problem: 
clearly the question is to be addressed by looking in detail at different models. 
In considering correlations between observables at a given epoch or in similar 
environments, the details of the radial modulus or other light scalar dynamics 
are not directly relevant. 

3.3 Scherk-Schwarz revisited 

As mentioned in the Introduction, Scherk-Schwarz symmetry breaking, in com- 
bination with varying radii, appears to lead to a relation between perturbative 
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gauge couplings and the mass scale of symmetry-breaking very different from 
radiative or nonperturbative sources of SSB. 

However, there are subtleties in the analysis which need to be addressed. 
It has recently been shown that Scherk-Schwarz breaking in a wide variety of 
situations can be reformulated as spontaneous breaking with v. e. v.'s being non- 
trivial functions of position in the extra dimension(s) (see e. g. |36jl. This may 
blur the distinction between Scherk-Schwarz breaking and breaking by radia- 
tive or nonperturbative means. Various connections have been made between 
Scherk-Schwarz supersymmetry-breaking in the fifth dimension and nonzero 
V. e. v.'s of auxiliary fields which might arise dynamically, both in heterotic 
M-theory simpler orbifold contexts [jg j . 

Such a connection might justify the choice of an extremely small dimension- 
less parameter a ~ 10^^'^ in to characterize the "twisting" of the boundary 
conditions. Thus what originally appeared as a geometric and explicit breaking 
induced by a nondynamical choice of boundary conditions, turns out to be a 
dynamical spontaneous breaking whose size may be determined by nonpertur- 
bative gauge theory effects. (In both cases the breaking is soft.) 

Hence, it is crucial to trace the mechanisms of symmetry-breaking back to 
their source, since consistent predictions for varying couplings depend on know- 
ing how the mass scales and v. e. v.'s of the model are generated. In the case just 
mentioned, since gaugino condensation would take place at a scale much lower 
than that of the inverse radius, we must take account of the variation in the 
condensate resulting from a variation of the radius. Thus the assumption that 
the S-S parameter a is constant, which would imply that the SUSY-breaking 
masses vary strictly with 1/R, is likely incorrect. However, if a S-S breaking 
parameter is a number of order 1 or quantized, we take it to be constant since 
it does not require a separate dynamical explanation. (See for a discussion 
of dynamics that could fix the value of such a parameter) . 

3.4 Units and definitions 

For varying extra-dimensional radii, the relevant dimensionless quantities are 
the radius relative to the cutoff of the extra-dimensional field theory i?A o , and 
relative to the extra-dimensional Planck length RM^\ Since by assumption 
(4 -|- (5)-dimensional physics above is invariant, we expect Ad/M^^ to be 
unchanged, hence the two ratios carry the same information. We will not be 
dealing with gravitational effects, hence we take RA^ as the "order parameter". 

In contrast to the usual designation of Mp as a constant, resulting in "Planck 
units" if we have Mp = 1 (other conventions import factors of ^/Stt, etc.), we 
could choose to work in "radius units" in which R is constant and Api and Mp 
vary. Thus "varying R" , in addition to affecting the classical relation between 
couplings in 4 -|- 5 dimensions and D = 4, implies RG evolution in the extra- 
dimensional theory, since we are changing the ratio of the UV scale to the IR 
scale(s). Since there are a large number of thresholds in such theories, with 
masses determined primarily by R, radius units may be more convenient. 
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3.5 RG "running" in 4 + S dimensions 

As is well known ,33 , the behaviour of coupling constants under change of 
scale is radically altered above the compactification scale if fields propagate 
in extra dimensions. The usual procedure is to decompose extra-dimensional 
operators into K-K modes whose masses are fixed multiples of and treat 
the massive modes as thresholds at which the effective four-dimensional beta 
function changes. With a large number of such modes, one finds an "averaged- 
out" behaviour of power-law running of the effective 4d coupling. Obviously, 
the running depends on which fields are allowed to propagate "in the bulk", 
and proceeds as normal in 4d if all fields coupled by a particular operator are 
"on the brane" . 

A cutoff is required since gauge and Yukawa couplings quickly become non- 
perturbative with increasing energy, being irrelevant operators in higher di- 
mensions For boundary conditions imposed at a scale Ad to be meaningful, 
the theory should be perturbative at this scale. By assumption, the coupling 
strengths of all operators in the (4 + 5)-dimensional theory at and above do 
not vary. 

We quote the relevant formula from j33j 

„,V-).„-'(A.,-|.n£4.„|^.gi|(^) -1) (17) 

where Ms = R^^ is the mass of the first Kaluza-Klein mode. Thus 

ai(/x-) y I y 27r 27r \Ad ) ) Ad 2t, Ad' 

(18) 

recalling that ai(Ao) oc {Ms/AdY', where the i'th gauge group propagates 
over 5i extra dimensions. The first term inside brackets arises from this classi- 
cal dependence, the second from logarithmic K-K contributions and the third 
from the power-law correction; the last term being the usual D = 4 scale de- 
pendence. Formally, the second and third terms inside brackets are suppressed 
by a loop factor relative to the first, however the power-law term may have a 
large enhancement, essentially due to a large number of K-K modes. However, 
we can choose the cutoff such that the classical term dominates, as follows. 

If couplings are perturbative at the cutoff Ad, we can choose a lower effective 
cutoff A', integrating out all modes between the two. For A' > Ms the form 
of Eq. H17() will remain the same and the extra-dimensional gauge coupling 
at A' can equally be taken constant. Thus we lower A' to just above Ms, 
making the power-law contribution of K-K modes as small as possible and, in 
practice, negligible. The quantum terms in the brackets can then be neglected as 
formally of higher order, since we do not expect b to be very large. We similarly 
neglect (the variation in) the difference between ln(;U~/A') and \n(^/Ms), and 

^^For energies far enough above , the behaviour of the compactified theory should 
reproduce that in uncompactified = 4 + 5 field theory, and the Kaluza-Klein modes can be 
approximated by a continuum. 
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set ai(A') ~ ai{Ms) in all static quantities. Whether {bi /In) A ln{n~ /Ms) is 
also of higher order, will depend on how the low energy scale is dynamically 
generated: if varied with the QCD scale Ac it could compete with the first 
term. 

Similar considerations for Yukawa couplings imply that the simple classical 
formulae 

A In a, (M^) = <5, A In ^ , A In y (M^) = + c) A In ^ 

are good approximations in the perturbative regime, where the Yukawa coupling 
y couples py = 0, 1, 2, 3 fields (either fermionic or scalar) propagating in extra 
dimensions and (3 — py) fields localized "on the brane" . The integer c takes the 
value for Yukawa couplings localized in extra dimensions, as in some models 
based on = 2 supersymmetry in D = 5, and 1 for a bulk coupling in the case 
where all fields propagate in the extra dimension. For a supersymmetric Higgs 
mass term nsHiH2 localized on the branc with Higgs in the bulk, us varies 
directly proportional to Ms, fig /Ms =constant; for other configurations (i. e. 
universal extra dimensions or localized Higgses) ns arises directly from a mass 
term in the D = 5 theory and thus varies as /jls/Ms oc {Ms/A')~^. 

We take one varying dimension, the generalization to more than one is 
trivial. Substituting for ai(Ms) and neglecting the variation in low-energy 
thresholds, including fi^ , relative to Ms, we find 

(19) 

Then using the naive formulae Aa ^ = Aay^ + Aag ^ and Aln(mp/me) = 
(27r/9a3(M5)) A\iia3(Ms), and neglecting running of a below the electroweak 
scale and possible variation in rUe/Ms, we find 

Alna = (5YCOs''ew + 52sm^9w + ^^^{5YbY + 52h2)ln^\A\n^^, 

\ ZTT Ms J A' 

Al 27r<53 / ^^ &3 , Mz\ . , Ms 

Aln^ = -9-(^«3 (^^) + 2^1-M^J^l" AT (^0) 
thus without superpartners we have 

* (^•'' - ^ W) '"^^ + * + (f - ^ 

(21) 

and with superpartners (in the case Ms 3> m) 

Rf^Ss (^5.9 - ^ In ^ j (^6y cos^ 0w + 62 sin^ dw + ^ {USy + ^2) In ^ 

'(22) 

Enforcing dy = 62 = S3 as required by an extra-dimensional GUT we obtain 
5.9-Iln(M./M,) 5.9-iln(M./M,) 
1 + 8D4 • T HMz/Ms) ^' 1 + • 12 HMz/Ms) 

(23) 
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If we use the more detailed treatment of the proton mass Eq. (|14|) then the 
value of R is simply reduced by a factor 0.78. 

In the small radius limit Ms — > Mq where the K-K modes do not influence 
unification we recover the SUSY-GUT expectation R ~ 33. In the opposite 
limit where Ms approaches Mz we would obtain ^ > 6. 

The suppression of R compared to its value in 4D GUTs was named k in j^l] , 
who found k ~ 0.16 in the limit of "low scale unification". Our estimate is 
At > 0.18 taking the value -R4DGUT ~ 33. Note that the estimate -R4DGUT ~ 36 
used in [HI] can also be obtained from our formula, after accounting for the 
running of agm below Mz through Alna = {a/a{Mz))Alna{Mz)- 

The above expressions for R are obtained without knowing any details of 
the supposed extra-dimensional GUT theory: the identity of the GUT group, 
the pattern of breaking, the localization of matter, etc. We do not even need 
to have a unified theory! The only parameter is the ratio Mz/Ms- Taking a 
value Ms = 5 TeV slightly larger than the minimum, without superpartners, we 
obtain R ~ 9.2: the logarithmic variation of a^^{Ms) with increasing ag ^(M^) 
is significant at low compactification scales. 

In more general "brane world" models, one has more freedom to choose the 
integers 5i and py. Clearly, the choice ^3 = leads to A In ^ = 0, up to effects 
(so far set to zero) of thresholds below Ms and variation in nhe/Ms- This is 
just the obvious scenario in which the strong force does not "feel" the variation 
because it is not propagating round the varying dimension. Data will turn 
out to favour small values of R, thus assigning iJs = is an obvious choice, if 
somewhat arbitrary in the absence of a concrete model. 

To assess the influence of thresholds and varying m^jMs we will now turn 
to more concrete models of mass generation (which, however, will not have the 
feature that ^3 = 0). 

3.6 Electroweak breaking by boundary conditions 

The possibility of electroweak symmetry-breaking triggered by Scherk-Schwarz 
supersymmetry-breaking on an orbifold j4()| I41j leads to a class of predictive 
scenarios where the variation of all observables, under a variation in radius, 
can be found with very few assumptions. The variations in gauge and Yukawa 
couplings (and /^s) as sketched above can be applied to obtain the behaviour 
of vu and fermion masses. The main result will be the variation of rrie/Ms 
and the resulting effect on /i. As we will see, given the assumption that all 
mass parameters in the electroweak sector originate from compacification, the 
fractional variation in fermion and vector boson masses relative to Ms cannot be 
more than a few times that in Ms/ A' (compared to around thirty times, in the 
scenarios of high-scale models with EWSB through dimensional transmutation) . 
Thus variations in charged thresholds are indeed small enough to neglect in 
these models. As for supersymmetric partners, the nature of SUSY-breaking 
ensures that the leading behaviour of soft masses is a constant times Ms, thus 
the variation of their thresholds is also negligible. Since (depending on details 

^^Of course, if unknown intermediate thresholds appear below Ms, all predictions are off. 
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of the model) superpartner masses are around the compactification scale, we 
use the SM beta-functions to run up to Ms (as in Eq. (^H)). 

As a first estimate, we truncate the Higgs potential to the quadratic and 
quartic terms and assume that first, the quadratic Higgs mass term is generated 
solely by the radiative effects of the top and stop loops and thus varies as 
M|y^, and second, that due to the underlying SUSY structure the quartic 
Higgs couplings vary as a2{Ms) + ai{Ms). 

Thus, in the one-Higgs model of in this approximation, the Higgs 
v.e.v., which varies as (|m|^|/A)^/^, behaves as /S.h\{yH /M^) ~ Aln(M5/A') 
and fermion masses vary as Aln(mg^i/M5) ~ (5/2)A ln(M5/A'). In the two- 
Higgs model of the Higgs v. e. v. is directly proportional to M5, A \d.{vh /Ms) ~ 
0, and the quark and lepton masses have a weak dependence lS.\D.{mq^i/Ms) ~ 
(1/2) A ln(M5/A'). The difference is simply due to different localization of Stan- 
dard Model fields. 

3.6.1 Detailed treatment of the Higgs potential 

One may treat the models more precisely by using the full Higgs potential, which 
may include terms besides the usual quadratic and quartic ones. In |41j the 
resulting dependence is = (const.) ((7^ + g''^)~^ R~^: the expected dependence 
on yt does not appear, due to a quirk of minimization of the potential, and 
the dependence on the quartic D-term coupling is different. Hence we have 
A\ii{vh/Ms) ~ (-l/2)Aln(M5/A') and /^\n{mq^i/Ms) ~ Aln(M5/A'). 

In the model of [101 the minimization conditions do not give a simple form 
for the dependence of vh on input parameters, due to the presence of two Higgs 
doublets and radiative corrections. The ratio of the two Higgs v. e. v.'s tan/3 is 
around 37 so we neglect 1 / tan^ /3 next to 1 and obtain 

A In ^ = A m I f V + f 2A, + 'I \ i'JL] ' \ (24) 

where is a quartic coupling proportional to yf with a weak logarithmic depen- 
dence on RMz cc Mz/Ms. Then clearly the dependence ns/Ms oc (M^/A')"^ 
is the dominant effect for the range of values of /ig considered in ^Uj. For 
Scherk-Schwarz parameter uj = \/2>, fis = 600 GeV and R^^ = 9TeV we find 
I\hi{vH/Ms) = 44Aln(M5/A') while for ^13 = 200 GeV and R'^ = 3TeV 
we find l^\n{vH/Ms) = 7Aln(M5/A'). Roughly, l^hi{vH/Ms) scales with 

The extreme sensitivity to /X5 and the fact that this parameter is put in 
by hand prevent the model from making any definite prediction. Essentially 
there is no dynamical explanation of the mu-term, it is a brane mass term in 
the D = 5 theory which we have to assume is non-varying relative to the extra- 
dimensional cutoff. It is not correlated with the compactification scale (times 
some coupling constants), thus, in radius units, it varies inversely compared to 
the radiative corrections triggering EWSB. If fis is somewhat larger than Mz 
it completely dominates the variation of vjj giving a large and arbitrary result. 

In general, if there are comparable mass scales with different dynamical 
origins, the variation of a quantity obtained by adding different mass terms can 
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be very sensitive to the proportion of different terms, thus the model parameters 
need to be known rather precisely. If the mu-term were to arise from the 
compactification scale then the dependence on the actual value of ns/Mz would 
be much milder and the scenario could be predictive. 

3.6.2 More precise estimate of R 

Finally we combine the variation of vh and fermion masses in the model of |41j 
with the detailed treatment of the proton mass to obtain the estimate of R. In 
this case the variations in nif/Mg is indeed small enough that one can neglect all 
threshold corrections and quark mass contributions to the variation of nip /Ms . 
Taking account of the variation Aln{me/Ms) ~ Aln(M5/A'), which simply 
subtracts 1 from the denominator of Eq. 1231 we find 



It has been noted that a divergent Fayet-Iliopoulos term might arise in models 
of this type the question of whether this could have a considerable effect 
on R is interesting but beyond the scope of this paper. 

3.7 GUT breaking by boundary conditions 

At the other end of the scale of compactification radii are scenarios where 
the energy scale of the extra dimension is just below the SUSY-GUT scale, 
as required if a GUT is to be broken by orbifold boundary conditions j441 
143) . These types of theory do not have dramatic low-energy predictions of 
detectable K-K modes, however they are claimed to be an improvement over 
four-dimensional GUTs in terms of predictivity and agreement with current 
values of sin^ Oyy and 03 (M^), at least for certain simple models. The ratio 
of the compactification scale to a fundamental cutoff may also play a role in 
generating fermion mass hierarchies, thus the variation of fermion masses and 
yt with a varying radius may be under control (although a full model explaining 
the ratio me/rrir has not appeared). 

However, as previously noted the mechanism of EWSB is not fully accounted 
for in such models: the dynamical origin of the hierarchy My/ /M^ is not speci- 
fied. Since this is also the major cause of uncertainty in four-dimensional GUT 
models, extra-dimensional GUTs of this type currently do not have predictions 
distinguishable from 4D GUTs for varying couplings. 

4 Comparison with data and conclusion 

There are four types of observable that may be considered in the redshift range 
of interest: a, X — > a^gprne/mp, Y — > a'^Qp and /u — > mp/me where — > signifies 
"is a measurement of". Current values are as follows: 



i? = 6.3±l. 



a-lAa = (-0.543 ±0.116) X 10"^ 0.2 < z < 3.7 [H 
X-^AX = (0.7 ± 1.1) X 10"^ z = 1.78 m 



(25) 
(26) 
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y^^Ay = (-0.16 ± 0.54) X 10"^ z = 0.68, (27) 
(-0.20 ± 0.44) X 10"^ z = 0.25 PJ (28) 
^u-^A/i = (5.0 ± 1.8) X 10~^ z = 3.02 (29) 

Theory and data can be usefully compared for a and /i; a variation in Qp cannot 
reliably be related to more fundamental parameters 

With three independent types of measurement and four observables, one 
may in principle check consistency and obtain values of gp or mp/nie in two 
different ways. Given measurements at similar redshifts and ignoring compli- 
cations due to spatial variation, one could eliminate gp as follows: 

Y Xfi 

^ ^ 5p, ^ ^ gp- (30) 

If the resulting fractional variation were consistent but large compared to a^^Aa, 
it would indicate a large quark mass contribution to gp, which would be rather 
unexpected. Similarly one can compare the "direct" measurement of /i to that 
extracted from 

y rup 
X rrie 

which will be consistent if and only if the values for gp are. The value of a itself 
derives only from the "direct" measurement. 

If we pretend that current data can legitimately be combined, then Y/a^ 
would give Alngp = (0.93 ± 0.6) x 10~^ whereas X^ijc? gives Aln^p = (6.8 ± 

2.1) X 10~^. The two values differ by over 2.5 standard deviations, so one 
hesitates to combine them; the second value also deviates from zero by over 3 
sigma. The variation in mp/me derived from Y/X would be A In// = (—0.9 it 

1.2) X 10^^, also differing from the "direct" measurement by over 2.5cj. But such 
combinations of measurements are strictly not correct, since they correspond 
to widely differing redshifts and environments. Such checks can only be taken 
seriously if the number and reliability of data increase significantly and the data 
sets truly overlap. 

However, both the measurement of /x at z = 3.02 and the derived value 
Y/X at redshift around 1 bound any model where smooth and not too rapid 
time variations in a imply variations in mp/me, since a nonzero variation in a 
appears to occur at least for redshifts 1 to 3. To stand a chance of agreeing with 
experiment, Aln/i should be in the range (—2.9,5.8) x 10~^, thus R should be 
in the range (—10.5, 5.5) (rounded to the nearest half integer). 

It is difficult for many of the models to fit in this window: the "best fit" 
among 4d unified theories discussed above is the case where (5^ and (is are 
around 34 (some scenarios discussed in may also survive). With some 
"adjustment" of (3^ to take account of a varying top Yukawa in concrete models 
of radiative EWSB a better fit might be obtained, but this hardly constitutes 
a motivation for further investigation. Our main result is that in models with 

^■^It has been claimed that cannot vary substantially, being essentially a Clebsch-Gordan 
coefficient flfci) . but also that it may have a substantial but unknown correlation with the 
strange quark mass |22|. 
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high-scale unification the variation of thresholds and fermion masses can be 
a leading contribution to the variations in a and fi: we presented a detailed 
discussion of these contributions. 

In extra-dimensional models, the variation of a compactified dimension rel- 
ative to the cutoff can lead to definite predictions, both for coupling constants 
and mass ratios, even if there are no unification relations. One only needs to 
know the physics below the compactification scale. Predictions depend mostly 
on the compactification scale Mg and on which fields propagate in extra dimen- 
sions. As a special case, we confirm a previous result [31] that the variation 
of the strong coupling is suppressed relative to the variation of a in extra- 
dimensional GUTs with low unification scale. 

If, as in some "brane-world" models, the gluons of QCD do not propagate 
round the varying dimension, while one or both of the SU(2)xU(l) groups do, 
the variation of the RG invariant scale Ac and the proton mass (relative to 
Ms) derive only from thresholds and quark masses and are likely to be small, 
while the variation in a is not suppressed. Thus one contribution to R is small. 
However, such a scenario also requires a concrete mechanism of EWSB to be 
self-consistent. If the variations of the electron and quark masses also turn out 
to be small, one would reach the region of R allowed by observation. 

Note that the coupling of a light scalar (p, supposed to be the source of 
the variation, to matter depends on the function dm/dcj), where m is the mass- 
energy of the body considered in Planck units. This quantity, which is cor- 
related with some of the deviations from GR, is at leading order given by 
d{mN /Mp) / d(f). Violations of the WEP due to composition-dependent forces 
will also include a term varying as {d/d(f)) (nip + nie — ran ) / Mp . Thus , scenarios 
in which the variations of rUp/Mp and vh /Mp are both small are better placed 
to satisfy bounds on both the variation of /i and on deviations from GR, than 
if R is small but nie/Mp and rUp/Mp both have large variations. 

In the extra-dimensional models we considered which do include a definite 
recipe for EWSB, all gauge bosons are in fact required to propagate round 
the compact dimension. In one model, the variation of vh relative to Ms was 
extremely sensitive to the undetermined value of the mu-term thus we were 
not able to arrive at a prediction. In another model with fewer parameters, 
agreement with data is tenuous in that the prediction for R is positive and 
larger than 5, even including the variation in nie/Ms. However it is possible 
that slightly different models of EWSB might result in definite predictions with 
a larger variation in me/ Ms which could bring R within an acceptable range. 
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